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ABSTRACT 

About one third of early-type barred galaxies host small-scale secondary bars. The formation and evolution 
of such double-barred galaxies remain far from being well understood. In order to understand the formation of 
such systems, we explore a large parameter space of isolated pure-disk simulations. We show that a dynamically 
cool inner disk embedded in a hotter outer disk can naturally generate a steady secondary bar while the outer 
disk forms a large-scale primary bar. The independent bar instabilities of inner and outer disks result in long- 
lived double-barred structures whose dynamical properties are comparable with observations. This formation 
scenario indicates that the secondary bar might form from the general bar instability, the same as the primary 
bar. Under some circumstances, the interaction of the bars and the disk leads to the two bars aligning or single, 
nuclear, bars only. Simulations that are cool enough of the center to experience clump instabilities may also 
generate steady double-barred galaxies. In this case, the secondary bars are “fast”, i.e., the bar length is close 
to the co-rotation radius. This is the first time that double-barred galaxies containing a fast secondary bar are 
reported. Previous orbit-based studies had suggested that fast secondary bars are not dynamically possible. 
Subject headings: galaxies: formation — galaxies: evolution — galaxies: structure — galaxies: stellar content 
— galaxies: kinematics and dynamics 


1. INTRODUCTION 

Double-barred (S2B) galaxies, consisting of a small-scale 
secondary bar embedded in its large-scale p rimary counter¬ 
part, w ere first described nearly 40 years ago (Ide VaucouleursI 
Il975h . Statistics of the fraction of S2Bs a mongst early-type 
galax i es have been obtained fro m optical (lErwin & SoarO 
I 2 OO 2 I: lErwinl[200^ and infrared (iLaine et al. 1 120021) observa- 
tions which showed that about one third of early-type barred 
galaxies are S2Bs. We still lack systematic surveys of later 
Hubble types because of t heir stronger dust extinction, espe¬ 
cially in the centr al region (lErwi ^ l2005h 
Observations dButa & Crocked Il99l: lEriedli & Martinetl 
ll993l:ICorsini et al.ll2003h have shown that the two bars rotate 
independ ently, which was expected from numerical simula- 
tions te.g. lShlosman & Helleil200^:lDebattista & Shenll2007L 
hereafter DS07). I n general, nested bars ca nnot rotate through 
each other rigidly (iLouis & Gerhardll988h . In the potential of 
two independently rotat ing bars, the orbits may not be closed 
in a ny reference frame. iMacieiewski & Sparkel (1 1 9971120001) 
and iMaciejewski & Athanassoulal (l2007h studied the orbits 
based on the concept of loops, which is a family of orbits 
whose population of particles return to the same curve, but 
not to the same position, when the two bars return to the same 
relative orientation. They also showed non-rigid rotation for 
loops. Dynamically decoupled secondary bars in S2Bs have 
been hypothesized to be a mechanism for driving gas past the 
inner Lindblad resonance (ILR) of the primary bars to feed 
the sup ermassive black holes that po wer active galactic nuclei 
(AGN) (iShlosman et 

Other works studied S2Bs from a purely kinematical point 
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of view. Ide Lorenzo-Caceres et al.l (120081) presented the 2- 
D stellar velocity and velocity dispersion maps of a sam¬ 
ple of four S2Bs, based on observations with the SAURON 
integral-field spectrograph. The high quality velocity disper¬ 
sion maps reveal two local minima, located near the ends of 
the secondary bar of each galaxy. They suggested that these 
cr-hollows appear because of the contrast between the veloc¬ 
ity dispersion of a hotter bulge and the secondary bar, as the 
secondary bar is dominated by ordered motion and thus has a 
low (J. 

The formation of S2Bs has been studied by numer¬ 
ical simulations. The best-known scenario for form- 
ing independently rotatin g double bars was proposed by 
lEriedli & Martinetl (Il993h : a pre-existing large-scale bar 
drives gas inflow into the central kiloparsec of a galaxy; 
Once sufficient gas has accumulated, it becomes bar- 
unstable and a dyn amically decoupled (gaseous) secondary 
bar forms (see alsolCombeslll994l^osman & Helled[2002: 
lEnglmaier & ShlosmanI 20041) . However, the S2B structures 
forming from gas are short-lived and gas dominated; this can¬ 
not explain the obse rved high abundance of S 2Bs in gas-poor 
early-type galaxies (iPetitpas & WilsonI 120041) . Ho wever, re¬ 
cently 3-D N-body-fhydrodynamical simulations bv lWozni^ 
(l2015h formed long-lived S2Bs that might provide a better de¬ 
scription of this process. The new stars form from central gas 
accumulations producing a dynamically cool inner disk. The 
stellar popul ations are in qualitative agreeme nt with the ob¬ 
servations of Ide Lorenzo-Caceres et al.r(l2013l) in that the sec¬ 
ondary bars are a few Gyr younger than their primary coun¬ 
terparts^ _ 

iRautiainen & Salol (119991) and iRautiainen et al.l (120021) re¬ 
ported the formation of S2Bs in purely collisionless studies, 
although these se condary bars often have a “vaguely spiral 
shape”. DS07 and lShen & Debattist^ (120091) performed well- 
resolved simulations of long-lived S2Bs with a pre-existing 
rapidly rotating pseudobulge without any gas. Einally, some 
simulations have indicated that the dark matter halo may 
sometimes play a role in generating S2Bs (iHeller et al.ll2007l: 
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ISaha & Maciejews]all2013h . Thus the conditions of S2B for¬ 
mation are still not well understood. Nevertheless, this va¬ 
riety of simulations shows that gas is not required to form 
secondary bars. 

Here we present new three-dimensional (3D) N-body sim¬ 
ulations that successfully generate S2Bs from simple and nat¬ 
ural initial conditions. We form long-lived S2B structures in 
pure-disk simulations by choosing different initial parameters 
for the inner and outer parts of the disk. Starting from differ¬ 
ent dynamical conditions, the inner and outer disks generate 
independent bar instabilities that result in a high probability of 
S2B formation. The paper is organised as follows. In Section 
[2|we describe the model setup. In Section [3] we summarize 
the results of exploring the parameter space. In Section [4] we 
show case studies of double-barred galaxies. A comparison 
with single-barred galaxies is presented in Section [5l Finally, 
in Section we summarize and discuss the implications of 
this work. 

2. MODEL SETUP 

The N-body simulations are e volved with a 3D cylin¬ 
drical polar gri d code, GALAXY (ISellwood & Vallurilll997l: 
ISellwoodI 120141) . All models consist of a live exponential 
disk (scale length stellar mass Mj) and a rigid dark matter 
halo, with no gas present. The simulations use units where 
G = Mo = Rd = V{) = To = \. The key point is that a dynam¬ 
ically cool inner disk is used to generate strong small-scale 
bar instabilities. As shown in Eq. [TJ the initial Toomre-Q 
profile is roughly constant at 2.0 in the outer disk (R> 1.75), 
while in the inner disk it decreases gradually to bq following 
a quadratic curve: 

{2.0-bQ){^)HbQ (7? <1.75) 

2.0 (7? >1.75). 

Thus the outer disk is dynamically hot, while the inner part 
of the disk is dominated by ordered motion. The only free 
parameter, Bq, can be used to set the dynamical temperature 
of the inner disk. The initial thickness zo is 0.1 all over the 
disk. 

In this study, the halo potential is logarithmic 

i>(r)=iy2ln> + r2), (2) 

with Vh = 0.6 and r/^ = 15. Because the models are very disk 
dominated, the halo density is low near the center where the 
secondary bar dominates. We expect that the angular mo¬ 
mentum transfer between the bar and the halo should be less 
important than mode ls in iDebattista & SellwoodI (120001) and 
lAthanassoulal (l2003h . Then rigid halos are useful for study¬ 
ing the complicated co-evolution of the two bars without the 
additional weaker evolution introduced by a live halo, by al¬ 
lowing high mass resolution in the nuclear regions. The ini¬ 
tial disk has 4x10^ equal-mass particles. A possible scaling 
to physical values is Mo = 4.0 x IO^^Mq and Rd = 3.0 kpc, 

which gives a velocity unit of Vb = km/s and 

a time unit of To = Rd/Vo — 12.3 Myr. The force resolution 
(softening) is 0.01 (corresponding to 30 pc). The forces in 
the radial direction are solved by direct convolution with the 
Greens function, while the vertical and azimuthal forces are 
obtained by fast Fourier transform. We use grids measuring 
Nr X N(f) X Nz = X 64 X 375 but have verified that increasing 


resolution does not affect our results by doubling Nr and A^. 
The vertical spacing of the grid planes was fe = 0.01. Time 
integration used a leapfrog integrator with a fixed time step 
5t = 0.04 corresponding to about 0.5 Myr. 

3. EXPLORING THE PARAMETER SPACE 



Fig. 1.— Exploring the parameter space of Bq and Mj. The long-lived 
double-barred, single-barred and nuclear-barred simulations are marked as 
“S2B”, “PB” and “NB”, respectively. The simulations whose two bars align 
with each other during the simulation are marked as “aligned S2B”. The black 
numbers at the right side of aligned S2Bs are the times when the alignment 
finishes. The shaded region shows the range of bq and when the inner 
disk suffers strong clump instabilities at the beginning. The red, green, and 
blue rectangles highlight the subsamples of standard S2B, clumpy S2B and 
NB models, while the dashed black rectangle marks the three models with 
the same disk mass to compare their evolution of their surface density. The 
results of the stochasticity test are shown as the fractions in red. 

Based on linear bar-formation theory (lToomrelll981h . the 
modes of bar formation are standing waves in a cavity, akin 
to the familiar modes of organ pipes and guitar strings. To 
form double-barred systems, we strive to build a disk with 
independent bar instabilities in the inner and outer regions. 
The responsiveness of the disk can be enhanced by decreas- 
ing Q, increasing surface density and reducing disk thickness 
(ISellwoodll 19891) Therefore, we explore the parameter space 
of Bq, Md and thickness. A set of simulations (not presented 
here) show that reducing the thickness of the inner disk does 
not strongly affect the formation of S2Bs. The results of sim¬ 
ulations with varying Bq and Md are shown in Figure [T] We 
classify the outcomes into five types based on their formation 
history and final morphology. There are three types of S2Bs: 
standard S2Bs, clumpy S2Bs and aligned (or coupled) S2Bs. 
Aligned S2Bs are unstable S2Bs whose two bars eventually 
couple into alignment after a few Gyr. The steady S2Bs form¬ 
ing from a violent clumpy phase when bQ< 0.45 are termed 
clumpy S2Bs. S2Bs not experiencing a clumpy phase are 
termed “standard” S2Bs. The two types of single-barred mod¬ 
els are divided into large-scale primary single-barred models 
(PBs) and nuclear-barred models (NBs) depending on the size 
of the bar. If the semi-major axis is less than 1.0, we classify 
it as a NB, otherwise, it is a PB. 

T he chaotic nature of disks leads to significa nt stochastic¬ 
ity (iMilled IT964 ISellwood & Debattistal 120091) . Because of 
the interaction of multiple non-axisymmetric components, the 
simulations presented here are more stochastic than typical 
simulations involving one bar. Therefore, we test the degree 
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of stochasticity for a subsample of models by changing only 
the random seed when generating the particle initial condi¬ 
tions. The results are shown as the fractions in red in Fig¬ 
ure [T] whose denominator corresponds to the total number of 
simulations we have run, while the numerator is the number 
of simulations forming steady S2Bs, including standard S2Bs 
and clumpy S2Bs. The fractions give an approximate proba¬ 
bility for forming long-lived S2Bs. When S2Bs fail to form, 
the outcomes can be either aligned S2Bs, or single bars, either 
NBs or PBs. 

4. THE FORMATION AND EVOLUTION OF DOUBLE-BARRED 
GALAXIES 

As shown in Figure[T] the most important condition for S2B 
formation from independent bar instabilities is given by the Q 
profile. When hq > 0.8, S2Bs do not form for a large range of 
disk mass. Instead most of such simulations form PBs only. 

4.1. A standard double-barred galaxy 

The model shown in Figure[2]is a standard S2B. The initial 
central dynamical temperature parameter, bq, is set to 0.5. In 
order to be sufficiently massive for bar instabilities, the mass 
of the disk is = 1.5Mo = 6.0 x lO^^M©. Figures [2] and [3] 
give an overview of the formation of the S2B over 100 time 
units (1.23 Gyr). The time evolution of the amplitudes of the 
primary bar (Aprim) and the secondary bar (Asec) is shown in 
Figured! 

From both analytical calculations and N-body simulations, 
we know that the length, strength and pattern speed reduction 
of bars are strongly affected b y the angular momentum ex¬ 
change with the out er disk fe.g. lLvnden-Bell & Kalnaj sill 9721: 

I Athanassoulal 120031) . Because of the low initial 2, the inner 
disk starts with strong m = 2 instabilities (r = 0 - 20, Figure 
[3]) with high pattern speed f^sec ^ 0.80. A small bar forms 
extending roughly to its co-rotation radius Rcr ^ 1. 0, which 
indica tes that it forms via the usual bar instability (iToomrel 
[ml. As time goes on, the small bar sheds angular mo¬ 
mentum and traps particles from the disk further out, which 
happens rapidly as the dynamical time scale is short in the 
inner disk. In this period {t = 0-50), Asec constantly in¬ 
creases, while the pattern speed of the secondary bar declines 
to (t = 40-60). 

The primary bar forms a little later and slower than the 
secondary bar through a quite similar process. It starts with 
^prim ^ 0.37 (t = 20-40) and declines to ~ 0.18 at r = 60-80. 
From the evolution of Asec (Figured)), we can see that the sec¬ 
ondary bar is significantly weakened by the formation of the 
primary bar, so that it can no longer extend to its co-rotation 
radius. It takes more than 300 time units for Asec to settle to 
a steady state, while the change of the primary bar amplitude 
is very small after t = 200. This S2B structure persists to the 
end of the simulation, lasting for more than 6 Gyr. The two 
bars are not rigid bodies when they rotate through each other. 
As seen in the insets of Figured! the dotted and dashed lines 
match the local maximum or minimum points of bar ampli¬ 
tudes well, and they also match the inflection points of Acj) 
curves in the inset of Figure [5! The secondary bar is stronger 
and rotates slower when the two bars are perpendicular, while 
it is weaker but rotates faster when they are parallel (see also 
DS07). Generally, the secondary bar amplitude and pattern 
speed variations are larger than in the primary bar. Such oscil¬ 
lating pattern spe eds and amplitudes are consis tent with loop- 
orbit predictions (iMacieiewski & Sparkell200Cih . 


In order to compare with observations, we use two meth¬ 
ods to measure the bar length. As shown in the top panel 
of Figure [5! the secondary bar peak amplitude of the stan¬ 
dard model is at R ^ 0.2. The length of the secondary 
bar, ^sec ^0.3, is measured by tracing half-way down the 
peak on the m = 2 amplitude plot, which is consistent with 
the value 0.4 given by the 10° deviation from a constant 
phase. The semi-major axis of the primary bar, a^nm ^ 3.0, 
gives a length ratio a^Qc/d^nm ^ 0.10-0.13 that is consistent 
with the typically obser ved length ratio of loca l S2B systems 
(median ratio ^ 0 .12) (lErwin & Sparkel [20Q1: lErwii]l2004 
iLisker et ani2006l) . The primary bar extends roughly to its co¬ 
rotation radius (7?c/?,prim ^ 3.0 ^ ^^prim), while the secondary 
bar is much shorter than its co-rotation radius (Rcr^sqc ^ 
1.5). If 1 .0 < Rcr I a S. E4, the bar is classified as a “fast” 
bar (e.g. iDebattista & SellwoodI l2000l: lAguerri et al.l 120031: 

ICorsini et~^ 120031 : IDebattista & WilliamsI l2004l)~ There¬ 

fore, the primary bar is a “fast” bar, while the secondary 
bar is “slow”. It is worth noting that the approximate co¬ 
rotation radius of the secondary bar (Rcr^sqc ^ 1-5) is close 
to the outer ILR radius of the primary bar (Rqilr, prim ^ 1-5), 
wh ich is considered as an evidence of CR-ILR coupling 
by iRautiainen & (119991) ( see also iPfenniger & NormanI 

ll990l:lFnedli & Martinetlll99^ 

Moreover, before the formation of the primary bar, the 
rapidly evolving inner disk exhibits recurrent transient spirals 
driven by the secondary bar as seen at ^ = 20 in EigurelD As 
the primary bar grows stronger, the nuclear spirals gradually 
disappear. Therefore, a reasonable inference is that the forma¬ 
tion of the primary bar efficiently suppresses nuclear spirals. 

4.2. A clumpy double-barred galaxy 

When bq < 0.45, all simulations experience a violent 
clumpy phase, shown by the shaded regions in Eigure [D As 
expected, with a colder inner disk, the clump instability be¬ 
comes more significant. The clump instability occurs mainly 
at R < 0.7 where the initial Q i s less than 0.8, whic h is roughly 
consistent with the criterion of iDekel et al.l (120091) {Q ^ 0.67) 
for clump growth in an unstable disk. A thinner inner disk 
will also strengthen the clump instabilities, but thickness is 
not as important as Q. When > 0.7, we never see clear 
clump instabilities even with a very thin inner disk. As shown 
in Eigure (T) the inner disk fragments into clumps in a short 
time. Because of clump-clump gravitational interactions and 
dynamical friction, the massive clumps sink to the center and 
merge to form a small-scale bar. The polar grid code we adopt 
concentrates spatial resolution at the geometric center of the 
grid. However, the formation of clumps results in the highest 
density not coinciding with the region of highest spatial reso¬ 
lution. Therefore, we have tested refining the grid by doubling 
the number of radial and azimuthal grid cells; the result of this 
test is similar to the result shown here. 

As with standard S2Bs, the two bars in clumpy S2Bs also 
form from independent bar instabilities of the inner and outer 
disks. Although the formation of the secondary bar is slightly 
delayed by the violent clumpy phase, as shown in Eigure [3 
the newly formed secondary bar can extend to its co-rotation 
radius as well, Rc/?,sec ^ 1-0. In our clumpy simulations, the 
clumpy phase always results in the formation of a small-scale 
bar at the beginning. However, they differ from standard S2Bs 
in terms of their dynamical properties. Eirstly, the amplitude 
evolution of the primary bar (Eigure O shows that it forms 
much later and more gradually in the clumpy S2B. Secondly, 
the iso-density contours (Eigure (T)) show that the secondary 
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Fig. 2.— The face-on images of the S2B simulation at various times, with isodensity contours superposed. The contours are equally separated in logarithm, 
showing the form ation process of the S2B structure over the interval 0 < r < 100. The surface density contours have been smoothed with an adaptive kernel 
JSilvermanll 19861) . The short and long straight lines mark the major axes of the secondary bar and the primary bar, respectively. A movie showing the evolution 
of the standard S2B in the primary bar’s corotating frame is available at hubble.shao.ac.cn/~dumin/S2BA^ideol.gif and the ApJ website. 



Fig. 3.— Power spectra of the m = 2 Fourier component as a function of radius. The signals of the m = 2 sectoral harmonic of the density are measured over a 
time interval of 20. The solid lines show fliR), the frequency of circular rotation. The dashed lines show Q + k/2 (upper) and Q-k/2 (lower), where k, is the 
epicyclic frequency. The two sets of solid and dashed lines correspond to the curves of the starting and ending times in each interval, respectively. Generally, the 
lower sets denote the starting times, while the upper sets denote the ending times. 
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Fig. 4.— Time evolution of the amplitudes of the secondary bar (top panel) 
and the primary bar (bottom panel) of the standard S2B. Apnm and Asec are 
dehned as the Fourier m = 2 amplitude over the radial ranges 0.96</?<3.01 
and R < 0.30, respectively. In the insets, the dashed lines mark the times 
when the two bars are aligned, while the dotted lines mark the times when 
they are perpendicular. 



Time 


Fig. 5.— Time evolution of the phase of the secondary bar, measured be¬ 
tween t = 510-560. The dashed straight line is the least-square ht, which 
gives the slope (T2sec}. The inset figure shows the phase difference, Ac/), 
between the phases of the bars and {Q)t, where (Q) is the pattern speed av¬ 
eraged over one relative rotation of the two bars; the thick line with a large- 
amplitude oscillation is for the secondary bar while the thin line with a small- 
amplitude oscillation is for the primary bar. The dashed (dotted) lines mark 
the times that the two bars are parallel (perpendicular) to each other as in 
Figure ID 



R 



R 

Fig. 6.— The m = 2 amplitude and phase as a function of radius at r = 500 
when all structures are stable. Here we take the simulations discussed in this 
paper as representative of each type of models to show their bar strength and 
length. The phase of the small-scale bars in the standard S2B, clumpy S2B 
and NB models are shown in the bottom panel. 

bar is more disky and central concentrated, while the sec¬ 
ondary bars of standard S2Bs tend to be more rectangular¬ 
shaped. We perform ellipse fitting for both models using 
I RAF, and find more negative B/[ for the secondary bar of the 
standard S2B. In Figure [TOl we show the central surface den¬ 
sity profiles of the three models. Because they have the same 
disk mass = 1.5, they share the same density distribution 
at f = 0. After the small-scale bar forms, at f = 80, the central 
surface density of the clumpy S2B is ^ 35% larger than both 
the standard S2B and NB model. From r = 80 to f = 500, the 
surface density of all these three models increases at the same 
rate. Therefore, the extra central surface density of the clumpy 
S2B is produced during the clumpy phase. Thirdly, the two 
bars of the clumpy S2B are both much longer than those of 
the standard S2B as shown in Figure [b] Measured by tracing 
half-way down the slope of m = 2 amplitude, the length ratio 
is a^Qc/a^nm ~0.8/5.0 = 0.16 that is slightly larger than stan¬ 
dard S2Bs. If the bar length is measured at the radius where 
the m = 2 phase deviates from a constant by 10°, is even 
larger at 1.0. Measured from Figure [8] (f = 200-220), Rcr, sec 
is about 1.3, then sec/^sec ^ 1.3-1.6. Thus the secondary 

bar is a quite “fast” bar extending close to its co-rotation ra¬ 
dius, which is quite common in our clumpy S2Bs. This sur¬ 
prising result demonstrates that it is dynamically possible to 
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Fig. 7.— The formation process of a double-barred galaxy from a clumpy phase. This simulation has = 1.5 and Bq = 0.3. 



Power spectra of the m = 2 Fourier component as a function of radius for the clumpy model in Figure [7] 


Fig. 8.— 

generate a stable “fast” double-barred system. As shown in 
the top panel of Figure [9l the change of Asec is quite small 
with the formation of the primary bar. Moreover, the two bars 
oscillate in the same way as standard S2Bs, but the oscilla¬ 
tion is smaller. These indicate the interaction of two bars is 
weaker. In Figure[7](f = 200,300), we can clearly see that the 
inner disk is dominated by the secondary bar. Thus the pri¬ 
mary bar is not efficient at trapping the secondary bar in the 
clumpy S2B models. A possible reason is that the torques 
from the primary bar are not strong enough to reduce the 
strength and length of the secondary bar, because the primary 
bar is much fatter and less massive than the ones of standard 
S2Bs. This is because the primary bar forms further out in the 
disk. Moreover, a more massive and concentrated secondary 
bar may also stabilize it against the primary bar. As shown in 
Figure [TOl in the very central region the clumpy S2B is more 
massive than the standard S2B. In conclusion, in our clumpy 
S2Bs, the primary bar is too weak to be efficient at trapping 
the secondary bar, in which case a “fast” double-barred sys¬ 
tem can be sustained for a long time. 

These clumpy S2B simulations demonstrate that a small- 


scale bar can be generated from a violent clumpy phase. Pre¬ 
vious N-body-fgas simulations showe d that clumps coalesce 
into a bulge who se shape was bar-like dMmegreen et al.ll2008L 
Figure 2 and 3). Ilnoue & Saitohl (120121) showed that clumpv- 
origin bulges are bar-like in their hydrodynamic (SPH) sim¬ 
ulations, and have exponential surface density profile, boxy 
(bar-like) shape and significant rotation. All these properties 
are consistent with the small-scale bars in our simulations. 
Therefore, these “bar-like” structures may be similar to our 
small-scale bars, which suggests that the secondary bar might 
often form from clump mergers. 

4.3. An aligned double-barred galaxy 

Figures [TT] and [T2| show the formation and coupling pro¬ 
cess of an aligned S2B. As shown in Figure [121 the initial bar 
instabilities form two bars with different pattern speeds be¬ 
tween r = 0 and 80. The oscillating properties are identical 
to those of the standard and clumpy S2Bs. During the period 
r = 80-270, the secondary bar gradually becomes shorter and 
weaker because of the interaction with the primary bar. At 
t ^ 280 (Figure [TT]) the secondary bar cannot pass through 
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Fig. 9.— Time evolution of the amplitude of the secondary bar (top panel) 
and the primary bar (bottom panel) of the clumpy S2B. Apnm and Asec are 
defined as the Fourier m = 2 amplitude over the radial ranges 2.00 </?< 4.07 
and R <0.96, respectively. 



Fig. 10.— The central surface density profiles of the standard S2B, clumpy 
S2B and NB models at = 1.5 highlighted by the dashed black rectangle 
on Figure [T] The black line is their initial surface density profile. 

the next perpendicular barrier where the rotational kinetic en¬ 
ergy of the secondary bar is minimal due to the gravitation 
torque from the primary bar. Then the pattern speed of the 
secondary bar decreases sharply to that of the primary bar. 
The secondary bar remains at fixed orientation relative to the 
primary bar during t = 280-300. The two bars share the same 
pattern speed, although they are still mis-aligned at a roughly 
constant angle. During this time, the surface density contours 
around the secondary bar become somewhat peanut-shaped. 
This shape transformation always happens in all our aligned 


S2Bs. Then the secondary bar falls back to alignment with 
the primary bar gradually, which is the lowest energy state. 
Sometimes, the secondary bar librates slightly around the ma¬ 
jor axis of the primary bar before complete alignment. Finally, 
they rotate at the same pattern speed until the end of the simu¬ 
lation. Generally, this coupling period occurs over a few tens 
of time units 100 Myr). 

Although the aligned S2B now resembles a single-barred 
galaxy, the relic of the secondary bar can be clearly seen 
as a nuclear peak in the m = 2 amplitude profile (Figure [Q). 
The peanut-shaped contours are another signature of aligned 
S2Bs that survive to the end of our simulations. Normally, the 
coupling process is dominated by the primary bar as the sec¬ 
ondary bar gradually becomes shorter and weaker. Only very 
few models show a coupling process dominated by the sec¬ 
ondary bar, in which the primary bar is nearly destroyed by 
the secondary bar during their coupling process. A fraction of 
the particles previously belonging to the primary bar are, in 
this case, trapped around the secondary bar, making it longer. 

For our standard S2Bs, the pattern speed ratio of the two 
bars, f^prim/^sec 5 varies from 0.36 to 0.55, while it is 0.25- 
0.43 for the clumpy S2Bs. Generally, the pattern speeds of 
the primary bars of the clumpy S2Bs are lower than of the 
standard S2Bs. Before alignment, the pattern speed of the 
aligned S2Bs is 0.35-0.49, i.e. similar to the range of the 
standard S2Bs. Therefore, it is hard to predict whether a S2B 
system will couple or not based only on its pattern speed ratio. 
Possibly, the mass and morphology of two bars also play a 
role. 

5. LARGE-SCALE SINGLE-BARRED GALAXIES AND 
NUCLEAR-BARRED GALAXIES 

When hq > 0.8, most of the simulations form large-scale 
single-barred galaxies, PBs, as can be seen in Figure [il Com¬ 
pared with the formation process of S2Bs, the bar instability 
starts with a lower pattern speed. By losing angular momen¬ 
tum to the outer disk, the bar gradually slows down. Then 
its co-rotation radius is pushed further out and more par- 
ticles at larger radiu s can be trapped forming a longer bar 
(lAthanassoulall2003h . In this way, they cannot form two in¬ 
dependently rotating bars. 

A fraction of simulations only leave small-scale bars after 
their weak primary bars are gradually dissolved, or, in some 
cases, never formed at all. These models are marked as “NB” 
(nuclear-barred model) in Figure [T] Although NBs are also 
a type of single-barred galaxy, compared with the PBs, they 
have different formation histories and properties. To form 
NBs, a mechanism is needed to prevent the mass of the outer 
disk from being trapped into the primary bar. The initial Q 
of the outer disk is about 2.0, which is close to the criterion 
(2.0-2.5) for non-axisymmetric stability in all disk mass dis¬ 
tributions (lAthanassoula & Sellwoodl[r986b . To quantify how 
sensitive the formation of the primary bar is to the change of 
2, we fix the Q profile of the inner disk, and vary Q in the 
outer disk from 2.0 to 2.3. \f Q> 2.2, the outer disk is un¬ 
responsive to bar instabilities, showing that bar formation is 
very sensitive to Q when Q is large. Even a small enhance¬ 
ment of Q in the hot outer disk can significantly suppress the 
formation of the primary bar. 

For the subsamples of NBs, standard S2Bs and clumpy 
S2Bs, we show the evolution of their azimuthally averaged 
2, (2), in Figure fT3l From t = 100 io t = 500, the changes 
in the inner disks are very small while the outer disks of 
the NBs become hotter than both standard and clumpy 
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Fig. 11.— The coupling phase of the two bars in the aligned S2B case. This simulation has Bq = 0.5 and = 1.1. Note the developing of peanut-shaped con¬ 
tours from t = 290. A movie showing the coupling behavior in the corotating frame of the primary bar is available at hubble.shao.ac.cn/~dumin/S2BA^ideo2.gif 
and the ApJ website 



Fig. 12.— Power spectra of the m = 2 Fourier component as a function of radius for the aligned S2B model in Figure FTTI 


S2Bs. Such dynamically hot outer disks become unsuitable 
to form a steady primary bar. In some cases, a weak and 
sometimes spiral-like primary bar forms in the early phase, 
but it is hard to stabilize for a long time. One possible 
explanation for the enhancement of Q is the spirals driven 
by the nuclear bar. A lot of studies have sho wn that spirals 
efficie nt ly heat disks to mor e eccen tri c orbits (jBertin & On 
1996 : iBarbanis & Woltieil 119671: ICarlberg & Sellwood 


1985: iBinnev & Lacev 119881: Jenkins & BinnevI 1199*5 
De Simone et all l2Q04 iMinchev & OuillenI I2Q06I) . In our 


simulations, we do see stronger spirals recur more frequently 
in NBs than in standard and clumpy S2Bs. 

6. DISCUSSION AND CONCLUSIONS 

We have been successful in forming double-barred galaxies 
with a set of simple initial conditions involving a cooler in¬ 
ner disk. By setting up an increasingly cooler disk towards 
the center, a S2B struct ure can form n aturally. Recent N- 
body-fgas simulations bv lWozni^ (12015h showed that the star 
formation in the central regions was responsible for stabiliz- 
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R 

Fig. 13.— The evolution of the mean value of Q of standard S2Bs, clumpy 
S2Bs, and NBs that correspond to the three groups of models in Figure [l] 
within the red, blue and green rectangles, respectively. The black lines are 
the initial Q profiles of these simulations. The error bars refer to the scatter 
of each group of simulations. 

ing the nuclear bars. New stars formed a dynamically cool in¬ 
ner disk whe re nuclear bar instabilities aros e, consistent with 
our re sults. iRautiainen & Salol (119991) and iRautiainen et al.l 
(I 2 OO 2 I) found that the secondary bar formed from the highly 
rotating central mode, and models with steeply rising rota¬ 
tion curves tended to have a long-lasting secondary bar. DS07 
used a pre-existing rapidly rotating pseudo-bulge to generate 
S2Bs. All these simulations suggest that a rapidly rotating nu¬ 
clear component might plausibly generate a small-scale bar 
with high pattern speed. The two bars form from indepen¬ 
dent bar instabilities of the inner and outer parts of disks. 
Thus the pa ttern speeds of the two bars d iffer significantly in 
such S2Bs. ISaha & MaciejewsO (12013h generated a S2B in 
a galaxy model whose gravitational potential was dominated 
by the dark matter halo. The two bars rotated at comparable 
pattern speed, but they developed from two independent bar 
instabilities as well. Therefore, we conclude that S2Bs can 
form by simply generating independent bar instabilities in the 
inner and outer disks, although the survival of S2B systems is 
stochastic because of the chaotic nature of disks. A dynam¬ 
ically cool inner disk embedded in a hotter outer disk is an 
easy and natural way to satisfy this condition. 

In the early phase of the formation of S2Bs, the secondary 
bars extend to their co-rotation radius, which indicates that 
they form from the usual bar instabilities as with the primary 
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bars. The angular momentum exchange among the hosting 
disk, secondary bar and primary bar determines the strength, 
morphology and final pattern speeds of the two bars. The pri¬ 
mary bars extend roughly to their co-rotation radius, so they 
are “fast” bars. For the standard S2Bs, although the secondary 
bars extend to the co-rotation radius at formation, they shrink 
to a much smaller size to become “slow” bars because of the 
suppression from the later-formed primary bars. The clumpy 
S2Bs show that the secondary bar can form from the violent 
clumpy phase. With a weaker primary bar, the clumpy S2B 
can host a “fast” secondary bar. Because of the interaction of 
the two bars and the disk, there are also a large fraction of un¬ 
stable S2Bs, with the two bars coupling to alignment in just a 
few Gyr. 

In our models, the secondary bars always form earlier than 
the primary ones, whi ch is consistent with pre vious collision¬ 
less simulations (^e.g. IRautiainen et ^120021 DS07). How¬ 
ever, because the whole galaxy is set up at the beginning, we 
cannot constrain the cold material accumulation of the cool 
inner disk. It is not entirely clear to us how a dynamically 
cool inner disk, required in forming S2Bs, may arise. It is 
generally thought that gas is driven from outside in (e.g. a 
possible channel is provided by the primary bar). In this case 
gas accumulated near the center can give rise to a cool disk 
that may be able to form a nuclear bar. However, an alterna¬ 
tive formation scenario is that the secondary bar forms from 
the violent clumpy phase in the early universe when the pri¬ 
mary bar is still not formed. In many N-body-fSPH simu¬ 
lations, these dumpy-origin small-scale bars could be easily 
mistaken for bulges. Therefore, the question remains whether 
the small-scale bar instabilities happens at the early time of 
the galaxy formation or after the formation of the primary bar. 
To better understand the formation of S2Bs, further numerical 
simulations are required. 
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